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We have carried out first-principles simulations, based on density functional theory, to obtain the 
atomic and electronic structure of (001) BaTiOs/KNbOs interfaces in an isolated slab geometry. We 
tried different types of structures including symmetric and asymmetric configurations, and variations 
in the thickness of the constituent materials. The spontaneous polarization of the layer- by- layer non 
neutral material (KNbOs) in these interfaces cancels out almost exactly the "built-in" polarization 
responsible for the electronic reconstruction. As a consequence, the so-called polar catastrophe is 
quenched and all the simulated interfaces are insulating. A model, based on the modern theory of 
polarization and basic electrostatics, allows an estimation of the critical thickness for the formation 
of the two-dimensional electron gas between 42 and 44 KNbOs unit cells. We also demonstrate 
the presence of an unexpected in-plane polarization in BaTiOs localized at the p-type Ti02/K0 
interface, even under in-plane compressive strains. We expect this in-plane polarization to remain 
hidden due to angular averaging during quantum fluctuations unless the symmetry is broken with 
small electric fields. 

PACS numbers: 73.20.-r,71.30.+h,77.80.-e,77.84.-s 



I. INTRODUCTION 

The surprising discovery by Ohtomo and Hwang^ of 
a metallic state at the interface between two good band 
insulating oxides, LaAlOa and SrTiOs, has triggered a 
large amount of new studies on polar oxide interfaces.^ 
Indeed, the quasi two-dimensional electron gas (2DEG) 
that forms when LaAlOa is grown on top of a Ti02 ter- 
minated (OOl)-surface of SrTiOa, i.e. when the inter- 
face between the two materials is LaO/Ti02, displays 
very different properties from those generated at inter- 
faces between standard III-V semiconductors (such as 
GaAs and A^Gai-a^As). Among them we find con- 
ducting carrier densities and electron effective masses or- 
ders of magnitude larger than those found at semicon- 
ductor interfaces.^ It is also fascinating how, depending 
on growth conditions, magnetic^ and superconducting^ 
ground states have been experimentally identified at this 
interface between non-magnetic insulating oxides. Very 
recently, two independent groups have proven how both 
magnetic and superconducting states might even coex- 
ist on the same sample,"^ a very unexpected result since 
magnetic order is usually considered detrimental to su- 
perconductivity. As a consequence of all these phenom- 
ena, interfaces in polar oxides can open the door to novel 
implementations of field effect transistors and to a new 
era of oxide electronics.!^ 

Despite this recent activity, many fundamental ques- 
tions regarding the origin and confinement of the 2DEG 
remain highly debated. Different models have been pro- 
posed to explain the experimental results. The pioneer- 
ing one invokes the so called "polar catastrophe"^ that 
arises from the polarization discontinuit}ff^ between the 
III-III polar LaAlOa film and the II-IV nonpolar SrTiOs 



layers along the [001] direction. Indeed, from the for- 
mal ionic charge point of view, LaAlOs can be described 
as a succession of positive (La^^O~^)^ and negative 
(Al+^02"^)"^ layers, while the alternating (Sr+^O-^)^ 

and (Ti^^O^^)^ layers of the perovskite structure of 
SrTiOs are charge neutral. But, aside this first ratio- 
nalization, other explanations can be found in the lit- 
erature for the origin of the 2DEG, among them: (i) 
the interlayer mixing between LaAlOa and SrTiOa and 
non-abruptness of the interfaces (with the formation of a 
few monolayers of metal lic La i-a^Sr^^TiOs^^), (ii) doping 
due to oxygen vacancie d^^ l ^^ * (including those produced 
in surface redox reaction^^^, and /or (i ii) the presence 
of charged defects and adsorbates.^^ All these mod- 
els highlight the importance of the growth conditions of 
these structures for the appearance and the behavior of 
the functional properties of the 2DEG. 

One of these properties, that is well reproduced by dif- 
ferent experimental groups on many samples grown with 
a variety of techniques, is the existence of a critical thick- 
ness, tc, in the number of layers of LaAlOa for the for- 
mation of the 2DEG. Thiel and coworkers^^ have demon- 
strated that, for the interfaces to be conducting, the num- 
ber of layers of LaAlOs has to be larger than four unit 
cells. The thickness of the polar layer increases up to five 
unit cells in n-type LaVOs /SrTiOs interfaces.!^ These 
observations are consistent with the fact that the conduc- 
tivity of SrTiOs-LaAlOa-SrTiOs heterostructures with 
dissimilar interfaces is reduced if their p- type (A102/SrO) 
and n-type (L aO/T i02) interfaces are spaced by less than 
six unit cells .^^^^ Remarkably, even below the critical 
thickness, a metal-insul ator t ransition can be driven by 
an external electric field.G^ These studies suggested the 
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possibility of the design of new polar interfaces where the 
appearance of the 2DEG could be switched on and off by 
the action of an external perturbation. 

Along this line, the replacement of one (or the two) 
materials at the interface by ferroelectric perovskites is 
particularly attractive. The spontaneous polarization 
present in these materials is very sensitive to electric 
fields and could be used to create a bound charge at 
the interface that could reinforce/deplete the 2DEG. Pre- 
vious theoretical works have been focused on I-V/II-IV 
interfaces [NaNbOs/SrTiOa,^^ and KNbOa/ATiOa in- 
terfaces, where A = Sr, Ba or PlJ ^^P^j . From the formal 
ionic charge point of view, I-V ferroelectric perovskite 
oxides, such as NaNbOa and KNbOa, are made of alter- 
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charged layers along the [001] direction (essentially as 
LaAlOs, although now the AO layers of the perovskite 
ABO3 structure are negative, while the BO2 layers are 
positive). Therefore, the layer- by-layer electrostatic of 
the previous I-V/II-IV interfaces is analogous to that in 
the LaAlOa/SrTiOa interface. Density functional the- 
ory simulations on non-stoichiometric (i.e. with a non- 
integer number of unit cells of the layer- by- layer non 
neutral perovskite), symmetric superlattices indeed sug- 
gested the existence of a 2DEG in KNbOs/ATiOa in- 
terfaces, switchable between two conducting states by 
the ferroelectric polarization orientation of the titanate 
layer.™ 

However, with the simulation boxes used in the previ- 
ous works only the n-type or p-type interfaces are present. 
It can be proved (see Sec. 4 of the Supplemental Mate- 
rial of Ref. 25) that, within this configuration, the local 
interface properties exactly reproduce those of the infi- 
nite isolated slab geometries, obviously beyond the crit- 
ical thickness for the formation of the 2DEG. In other 
words, the calculations of Refs. 23 and Ref. [24l show the 
charge distribution and properties after the electronic re- 
construction has taken place, but nothing is said about 
the magnitude of the critical thickness for the formation 
of the 2DEG.26 

In this work we carry out first principles calculations of 
BaTiOa /KNbOs interfaces where we explicitly avoid the 
issue of lack of stoichiometry in the simulation box. We 
have found that, due to the large lattice screening pro- 
vided by the KNbOs layer, the critical thickness for the 
formation of the 2DEG is one order of magnitude larger 
than in LaAlOa/SrTiOs interfaces. An unexpected re- 
sult of our simulations is that an in-plane polarization 
develops on the BaTiOa side of a (Ti02/K0) p-type 
BaTiOa/KNbOa interface even when BaTiOa is subject 
to in-plane compressive strains. We explain this effect 
using basic electrostatic arguments. 

The rest of the paper is organized as follows. After 
summarizing the basic theory behind the electronic re- 
construction in Sec. [llj we present the computational 
details used in our simulations in Sec. IIIII The first- 
principles results, together with the relevant comparisons 



to the model, can be found in Sec. IV 



II. BACKGROUND ON THE "POLAR 
DISCONTINUITY" MODEL 

In order to establish the nomenclature and the basic 
theory that will be used later, we review the most im- 
portant points of the "polar discontinuity" model. Al- 
though this model has been invoked since the discovery 
of the 2DEG at polar oxide interfaces,^ only recently it 
has been rigorously rationalized with explanations firmly 
rooted on the modern theory of polarization (for a re- 
cent review, see Ref. 27 and references therein). This has 
been developed by Stengel and Vanderbilt in Ref. [10] for 
insulating interfaces, and later generalized by Stengel for 
the case of a non-zero surface density of "free" charge 
in Ref. EH and to the case of surfaces in Ref. |25] The 
theory presented in these works is absolutely general, and 
we strongly point the interested reader to those milestone 
papers. Here, we particularize it to the conditions con- 
sidered in this work, and estimate the critical thickness 
for the formation of 2DEG in the case where any of the 
two materials forming the interface is ferroelectric. 

The standard nomenclature used in the literature of 
polar oxide interfaces denote the material that is non- 
neutral layer- by-layer (i.e. LaAlOs) as polar^ and the 
material that is neutral layer- by-layer (i.e. SrTiOa) as 
non-polar. Rigorously, this notation does not apply here 
since the two materials that constitute our interfaces are 
ferroelectric and, therefore, might undergo polar phase 
transitions with the appearance of a non- vanishing spon- 
taneous polarization. Nevertheless, for the sake of con- 
sistency with previous works we will maintain the con- 
vention and refer to the ferroelectric non-neutral layer 
by layer material (i.e. KNbOs) as the polar material 
and the ferroelectric neutral layer-by-layer material (i.e. 
BaTiOs) as the non-polar one. 

During the development of the model we assume a n- 
type interface, simulated within an isolated slab geome- 
try, with the non-polar material at the left and the polar 
material (with a formal ionic charge of ±e alternating 
from layer to layer, where e is the magnitude of the elec- 
tronic charge) at the right [see Fig.[T]^a)]. The generaliza- 
tion for other configurations is straightforward, changing 
the appropriate signs when required. 

Within the modern theory of polarization, we can com- 
pute the "formal" bulk polarization from the positions 
of the atomic nuclei and the center of localized Wannier 
functions. This decomposition of the charge (nuclear and 
electronic) into localized contributions allows for a simple 
classical interpretation of the bulk polarization in terms 
of a point charge model, and rescue the Clausius-Mossotti 
formulation. 

In the perfectly ideal structure without rumpling, 
where an atomically sharp junction in the absence of de- 
fects is supposed, all the atoms at a given layer lie at 
the same plane [Fig.jlja)]. Then, we can always choose 
unit cells that tile the crystal under appropriate primitive 
translations, and that leave the left-over interface region 
charge neutral.'^ (It is important to note that, for the 
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moment, we are assuming that the thickness of the po- 
lar layer is below the critical thickness for the formation 
of the 2DEG.) Then, the magnitude of the dipole of an 
individual (A0)-(B02) unit in this material is d = ea/2, 
where a is the out-of-plane lattice constant. The sign 
of the dipole is always directed from the negative to the 
positive layer, so in the considered configuration the sign 
is negative [see Fig. [TJa)]. This dipole corresponds to a 
"built-in formal" polarization (calculated for the previ- 
ous choice of unit cell for the primitive basis of atoms 
and Wannier functions) of 
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where Q is the volume of the unit cell of the polar material 
and S is the cell surface. Analogously, since the layers 
are formally charged neutral in the non-polar material, 
= n 

non— polar 

Now, we can wonder what would happen if the mate- 
rials that constitute the interface are ferroelectric, with a 
non-vanishing ferroelectric contribution to the polariza- 
tion (this might be also the case when a compressive 
epitaxial strain is applied to the SrTiOs/LaAlOa inter- 
face^. In this situation, the ferroelectric polarization 
must be added to the "built-in" polarization of the polar 
layer 

The difference in the polarization at the interface be- 
tween the two materials produces a surface density of 
bound charge. 
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Classical electromagnetic theory (Gauss's theorem) 
teaches us that this sheet of interfacial charge gives rise 
to a change in the macroscopic electric field in the two 
materials of 
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where eo is the dielectric permittivity of vacuum 
[Fig.jljb)]. The exact magnitude of the fields depends on 
the electrostatic boundary conditions, extremely linked 
with the geometry of the simulation boxes used in the 
computations (see Supplemental Material of Ref. [28] for 
a complete review). 

The field might induce strong structural changes in 
the materials and polarizes the electronic Wannier func- 
tions. Both facts translate into the development of a 
field-induced polarization AP [Fig. [ijc)] that tends to 
screen the discontinuity of the total polarization P, 
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and, consequently, of the macroscopic electric field. At 
the end, a self-consistent solution of Eqs. ([2|-(|4| is 
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FIG. 1. (Color online) (a) Schematic representation of a 
pristine polar interface in an isolated slab geometry. Atoms 
are represented by balls: O (small size, red), Ba (large size, 
yellow), Ti (medium size, blue), K (large size, green), and 
Nb (medium size, gray), with the corresponding atomic layer 
printed on top. The dots at the left of the BaTiOs indicates 
that we assume a thick layer. At the other side, a free KO- 
terminated surface of the polar KNbOs is assumed. Numbers 
below each layer indicate the formal ionic charge. The in- 
terface is marked with a red dashed line. Our choice for the 
unit cells that tile the entire crystal are represented by black 
dashed boxes, a stands for the out-of-plane lattice constant of 
the polar material unit cell, (b) Electrostatic potential of the 
polar interface in the isolated slab geometry. Within our elec- 
trostatic boundary conditions, the macroscopic electric field 
in BaTiOs is forced to be zero. Note that the field within 
the non-neutral layers points away from the interface, (c) In- 
duced polarizations AP due to the macroscopic electric fields, 
(d) Schematic representation of the energy bands of the polar 
interface. 



achieved where the polarizations in the two materials 
are in equilibrium with the corresponding macroscopic 
electric fields, and the total energy of the system is min- 
imized. 

For a sufficiently thin polar layer thickness, before the 
electronic reconstruction takes place, the absence of free 
charge at the interface requires the normal component of 
the electric displacement field D to be preserved. 
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The finite electric displacement is an input parameter of 
the model. In a first-principles simulations its value can 
be set by hand using the virtual crystal approximation 
to introduce external fractional charges in the surface 
atoms layers, while constraining the macroscopic electric 
field to be strictly zero in the vacuum region. Other 
authors fix the atomic positions on the surface unit cell 
to some specific values.*^ 

From the definition of the electric displacement field, 



-^polar — ^O^polar H~ -Ppolar 

= eo^polar + ^polar + ^polar + ^^polar, (6) 

and assuming that the polar material behaves as a linear 
dielectric of susceptibility Xpoiar around the spontaneous 
polarization structure, then 

^Ppolar ^0 Xpoiar Xpoiar 7 (7) 

and Eq. ([6| transforms into 

Xpoiar = ^0 Xpoiar Xpoiar + Xpoiar ^ Xpoiar' (^) 

where epoiar = 1 + Xpoiar is the dielectric constant of the 
polar material. From Eqs. ([5| and ([8|, 
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This electric field tilts the electronic bands of the po- 
lar layer [Fig. [ijd)]. At a given critical thickness, tc, the 
top of the valence band of the polar material reaches the 
level of the bottom of the conduction bands Beyond 
tc, a Zener breakdown takes place, with the concomitant 
transfer of charge from the surface of the polar mate- 
rial to the interface. The magnitude of tc can be easily 
computed from Eq. (|9| as 
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where A is the interfacial potential step. A will depend 
on the type of band alignment and on the particular in- 
terface (p or n). Here, according to the band alignment 
of Fig. [l[d) for the n interface is given h^p^ 



rppola: 
gap 



(11) 



where E^^^^ is the band-gap of the polar material. Note 
that the band alignment in our interfaces, where the band 



gaps of BaTiOs and KNbOa are very similar, is of type II, 
different from that in the prototypical LaAlOa/SrTiOs 
case (type I). For schematic views of the types of inter- 
faces according to the band offset see Ref. |321 



The first conclusion that can be drawn from Eq. (10) 



is that the polarizability of the polar layer is essential 
to determine the critical thickness for the formation of 
the 2DEG (it is directly proportional to Cpoiar)- The 
role played by the polar distortions to avoid the polar 
catastrophe in LaAlOa/SrTiOs interfaces has been con- 
firmed by first-principles simulations.^^ The importance 
of the extra screening due to lattice relaxations has also 
been discussed in the strongly related LaTiOa/SrTiOs 
inter faces. I^^M^ Indeed, the success for explain ing the 
critical thickness for the formation of the 2DEG,l^ii^^ 
together with the electrostrictive effect on the polar 
LaAlOa films,'^are between the most important achieve- 
ments of the polar discontinuity model. 

The second conclusion is related with the relationship 
between D and the ferroelectric polarization in the non- 
polar layer. In many works, the simulations try to re- 
produce the behavior of a thin polar layer on top of 
a thick non-polar substrate. In such cases, the macro- 
scopic field in the non-polar materials is forced to be 
zero, either by symmetry or by using a dipole correction 
in vacuum. As a consequence, D = Pnon-poiar- If the 
ferroelectric contribution to the polarization in the non- 
polar layer points in the same direction as the "built-in" 
polarization in the polar one, then it contributes to in- 
crease the critical thickness. This has been proven in 
Ref. l29l for the case of a SrTiOs/LaAlOs interface sub- 
ject to epitaxial strain [Eq. (10) is equivalent to Eq. (4) 
in Ref. 29 taking into account that, in this particular 
system P^^,^^ = P^^^Os = and I) = Pg^.Os]' 

Finally, the third conclusion is that there is also a 
strong influence of an eventual spontaneous polarization 
of the polar material in the value of tc. In particular, 
if Ppoiar close in magnitude to Ppoi^r points in 
the opposite direction, so the "built-in" polarization can 
be almost compensated by the ferroelectric contribution 
to the polarization, a large cancellation of the term in 
parentheses in the denominator of Eq. (10) is produced. 



with the concomitant increase in the critical thickness. 
In a previous work, Murray and Vanderbilt^^ have esti- 
mated how in SrTiOs/KNbOs superlattices the system 
would not become metallic until the number of layers of 
KNbOa is larger than 32. 



To further validate Eq. ( 10 ) in an isolated slab geome- 



try, we have carried out simulations on BaTiOs/KNbOs 
interfaces. The motivation for this choice is four fold: (i) 
under appropriate compressive in-plane strains, KNbOa 
is a ferroelectric polar material, with the spontaneous 
polarization pointing along the [001] direction, and with 
a "built-in" polarization of Xpoiar — ^KNbOs ~ e/2S ~ 
53 /iC/cm^ (value computed at the theoretical in-plane 
lattice constant of an hypothetical SrTiOs substrate, ay 
= 3.874 A), (ii) under this mechanical boundary condi- 
tion, both KNbOs and BaTiOa can be stabilized with 
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the same tetragonal PAmm symmetry, (iii) the theoret- 
ical spontaneous polarization of our tetragonal KNbOa 
in bulk is Ppoi^r = ^KNbOa = /iC/cm^, close to the 
built-in polarization, and (iv) we can also test to which 
extent the ferroelectric contribution to the polarization 
in the BaTiOa layer is dominated by the imposed value of 
the displacement field in the simulations. First-principles 
results and comparison with the previous model will be 



presented in Sec. IV 



III. COMPUTATIONAL DETAILS 

We have carried out density functional first-principles 
simulations based on a numerical atomic orbital method 
as implemented in the Siesta code.--^ All the calcu- 
lations have been carried out within the local density 
approximation (LDA), using the Perdew and Zungep^ 
parametrization of the Ceperley and Alder functional^ 
to simulate the electronic exchange and correlation. This 
choice avoids the systematic overestimation of the ferro- 
electric charactei^^ of perovskite oxides found in other 
commonly used functionalJ^H based on the generalized 
gradient approximation. This is an important point in 
this study, since the dielectric properties of the oxides at 
the bulk level might determine the behavior of the inter- 
faces, in particular a tendency for "overscreening" of the 
2DEG when the ferroelectric properties are favored. 

Core electrons were replaced by ab-initio norm con- 
serving pseudopotentials, generated using the Troullier- 
Martins scheme, in the Kleinman-Bylander fully non- 
local separable representation.^ Due to the large overlap 
between the semicore and valence states, the semicore 3s 
and 3p electrons of Ti, ?>s and 3p electrons of K, 45 and 
Ap electrons of Nb, and 5s and bp electrons of Ba were 
considered as valence electrons and explicitly included in 
the simulations. K, Ti, Nb and Ba pseudopotentials were 
generated scalar relativistically. The reference configura- 
tion and cutoff radii for each angular momentum shell 
for the pseudopotentials used in this work can be found 
in Ref. |46| for Ba, Ti, and O, and in Table [l| for K and 
Nb. 

The one-electron Kohn-Sham eigenstates were ex- 
panded iiia basis of strictly localized numerical atomic 
orbitals.ElHlweused a single-^ basis set for the semicore 
states of K, Ti, Nb, and Ba and double-^ plus polariza- 
tion for the valence states of all the atoms. For K (Ba), 
an extra shell of ?>d (bd) orbitals was added. All the pa- 
rameters that define the shape and range of the basis 
functions were obtained by a variational optimization of 
the eneigy^^ in bulk cubic BaTiOa (for Ba, Ti, and O), 
and of the entalph}/^ (with a pressure P = 0.03 GPa) in 
bulk cubic KNbOa (for K and Nb, the basis set of O was 
frozen to that obtained in BaTiOa). 

The electronic density, Hartree, and exchange correla- 
tion potentials, as well as the corresponding matrix el- 
ements between the basis orbitals, were calculated in a 
uniform real space grid. An equivalent plane wave cut- 



off of 1200 Ry was used to represent the charge density. 
For the Brillouin zone integrations we use a Monkhorst- 
Pack sampling^ ^ equivalent tol2xl2xl2ina five atom 
perovskite unit cell. 

To avoid the problem of artificially charging the in- 
terface with the use of non-stoichiometric superlattices 
with periodic boundary conditions, we follow the pro- 
posal of Lee and Demkov.^^ Within this approach, 
the calculations were performed on vacuum-terminated 
(KNb03)m/(BaTi03)^/(KNb03)m slabs, where the 
number of KNbOa cells, m, is always an integer number 
to guarantee the electrical neutrality of the interface. In 
particular, we have focused on two kinds of systems: sym- 
metric slabs where both interfaces between BaTiOa and 
KNbOa are of the same kind (either p-type Ti02/K0, 
or n-type BaO/Nb02 interfaces), and asymmetric in- 
terfaces where one interface is Ti02/K0 and the other 
BaO/Nb02. In both cases we have relaxed structures 
containing a different number of KNbOs and BaTiOa 
unit cells, given by the subscripts m and /, respectively. 
In the asymmetric interface a dipole slab correction is 
used to guarantee that the electric field in vacuum van- 
ishes. 

The in-plane lattice constant was fixed to the theoret- 
ical one of an hypothetical SrTiOs substrate (ay = 3.874 
A). Under this constraint, both BaTiOa and KNbOa are 
under conditions of compressive epitaxial strains. This 
will be relevant for the discussion of Sec. IIVBI 

Starting from ideal reference structures, built by piling 
up the corresponding unit cells of bulk strained materi- 
als without rumpling, the atomic coordinates are relaxed 
until the maximum component of the force on any atom 
was smaller than 0.01 eV/A. In the symmetric slabs, the 
minimization is performed in a two step process: First, 
mirror symmetry planes are imposed on the central layer 
of BaTiOs to avoid the development of a polarization in 
any direction. Then, symmetry is broken displacing co- 
herently the cations by hand, and a second relaxation is 
carried out without any imposed symmetry. In the asym- 
metric slabs, the relaxations are carried out in a single 
step, since the symmetry is broken directly in the initial 
reference structure. 

To stablish the notation, we will call the plane parallel 
to the interface the (x, y) plane, whereas the perpendic- 
ular direction will be referred to as the z-axis. 



IV. RESULTS 
A. Out-of-plane lattice relaxations and screening 

In order to characterize the atomic displacements 
induced by the relaxation, we define the "out-of- 
plane" rumpling parameter along z of layer i diS r]^ = 
[z{Mi) — z{Oi)] /2, where z{Mi) and z{Oi) are, respec- 
tively, the z coordinates of the cations and the oxygens 
at a given layer i [Fig.|2|a)]. We also define the "in-plane" 
rumplings {rj^^rf)^ in an equivalent way, as represented 



6 



TABLE I. Reference configuration and cutoff radii of the pseudopotential used in our study. Because of the inclusion of 
the semicore states in valence within the TrouUier-Martin scheme, K and Nb pseudopotentials must be generated for ionic 
configurations (ionic charge of +1). However, these are more suitable than the neutral ones, given the oxidation numbers of 
these atoms in the perovskites. Units in Bohr. 
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FIG. 2. (Color online) Schematic representation about the 
procedure to calculate (a) the out-of-plane rumpling, 77^, and 
(b) the in-plane rumplings, 7^^ ^if . Solid lines represent the 
position of the atoms in the ideal unrelaxed structure. After 
relaxation, the atoms move in the directions indicated by the 
arrows. Meaning of the balls as in Fig. [l] 



In all the studied systems a large out-of-plane rumpling 
is observed within the KNbOa layers, with a magnitude 
that is essentially independent of their thickness and the 
kind of interface: symmetric p [Fig. (sja)], symmetric n 
[Fig. |3jb)] or asymmetric [Fig. [sj^c)]. Two oxide layers 
away from the interface, the layer-by- layer rumpling con- 
verges to a rather uniform sawtooth pattern, with values 
slightly larger than those observed in bulk KNbOa un- 
der the same epitaxial conditions. This fact is consis- 
tent with liKNbOal ^ I^KNbOal' ^ith APKNbOs tending 
to compensate for the difference in order to screen the 
polarization discontinuity at the interface. Only in the 
neighborhood of the free surface, a small deviation from 
this trend is obtained due to the larger relaxations on the 
surface atoms. 

The dipole slab correction ensures that the displace- 
ment field in vacuum vanishes, = 0. Due to the ab- 
sence of surface external charges in the simulations, this 
value is preserved at the interfaces. This implies that a 
polarization in the BaTiOa layer induces a depolarizing 
field that is responsible for a large electrostatic energy. 
Therefore, under this electrical boundary conditions, no 



polarization is expected on BaTiOs. Both facts are well 
reproduced in our simulations, where we observe that 
the out-of-plane polarization vanishes within the BaTiOa 
layer in all cases. This happens even when the mirror 
symmetries are lifted by hand (in the symmetric inter- 
faces), or spontaneously (in the asymmetric slabs). 

According to the model developed in Sec. [TT| the po- 
larization in the polar layer tends to screen the disconti- 
nuity of the polarization at the interface and, therefore, 
its sign opposes that of the "formal" polarization. The 
lattice screening avoids the development of an electric 
field in the polar material that would result in the tilt- 
ing of its bands and, for sufficiently large thicknesses, to 
a Zener breakdown and accumulation of charge at the 
interfaces. The lattice screening provided by KNbOa is 
much stronger than the one anticipated in LaAlOs, where 
the structural distortion sustains the insulating behavior 
up to only 5 overlayers of LaAlOs on SrTiOs.'^ As a 
consequence, in our simulations all the computed struc- 
tures are insulating. The reason behind this is that the 
LaAlOa is a wide band gap insulator with a low dielec- 
tric constant (e^ = 25) and no ferroelectric instability 
(^LaAlOs ~^)' -'-t costs somc energy to polarize it. On 
the contrary, KNbOs is a ferroelectric oxide that polar- 
izes spontaneously and contributes to reduce the field 
[making smaller the numerator of Eq. (|9|], and increase 
the critical thickness for the formation of the 2DEG [mak- 
ing larger the denominator of Eq. (10)]. The reduction of 
the internal field within KNbOs can be directly checked 
from the nanosmoothed^^^ electrostatic potential in the 
slabs. Independently of the geometry, the magnitude of 
the field amounts to 0.024 eV/A [see red dashed lines in 
Fig. [4] to be compared with the roughly constant electric 
field of 0.24 eV/A ( one order of magnitude larger) found 
in LaAlOs/SrTiOs interfaces .^^^^ This is in very good 
agreement with the prediction of the electrostatic model 
developed in Sec. |ll| using a dielectric constant of bulk 
KNbOs around the ferroelectric structure under the same 
mechanical boundary condition of 25.0,^^ Eq. ^ yields 
a value of 0.023 eV/A. 

Since the interfacial potential steps between KNbOs 
and BaTiOs are 1.08 eV for the Ti02/K0 interface and 
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FIG. 3. (Color online) Out-of-plane (solid blue) and in- 
plane (long-dashed red) lattice polarization calculated for 
(a) symmetric-p [/ fixed to 4.5 and m — 2 (diamonds), 
m = 3 (squares), m = 4 (circles), m = 5 (crosses)], (b) 
symmetric-n (m=5, / — 4.5), and (c) asymmetric (m=4, /=8) 
BaTiOs/KNbOs slabs. Short-dashed (dot-dashed) lines rep- 
resent the rumplings of the Nb02 (KO) layers in bulk KNbOs 
under the same epitaxial constraint. The arrows point along 
the direction of the "built-in" polarization. Vertical dashed 
lines indicate the position of the interfaces. The in-plane rum- 
plings plotted here correspond to rjxy = +77^ defined in 
Fig.H 



1.31 eV for the BaO/Nb02 interface (both values com- 
puted using the recipe given in Ref. |55] when the bands 
of one material are tilted), the estimated critical thick- 
ness to trigger the polar catastrophe under the condi- 
tion of a vanishing electric displacement according to 
Eq. (10) are, respectively, 170 A and 180 A (between 



42 and 44 unit cells), one order of magnitude larger than 
in LaAlOa/SrTiOs. 



B. In-plane polarization 

Even though the ground state of bulk KNbOs and 
BaTiOs at zero temperature is rhombohedral, where 
the polarization displays both in-plane {x^y) and out- 
of-plane {z) components, when a compressive in-plane 
strain is applied the polarization in the (x, y)-directions 
is strongly reduced or even suppressed. In particular, 
when KNbOs or BaTiOs thin films are grown on top 
of a SrTiOs substrate, the in-plane polarization of these 
materials vanishes and the tetragonal c-phase is stabi- 




FIG. 4. (Color online) Macroscopically averaged internal 
electrostatic potential (black solid line) for (a) symmetric-p 
(m = 5, Z = 4.5), (b) symmetric-n (m = 5, / = 4.5), and 
(c) asymmetric (m=4, /=8) BaTiOs/KNbOa slabs. The red 
dashed lines indicates the region used to extract numerically 
the macroscopic field. 



lized.'^^'^ However, contrary to current thought, in our 
calculations of the BaTiOs/KNbOs interfaces, we can 
observe in Fig. [3] the appearance of a moderate in-plane 
polarization in these systems along the [110] direction. 
The layer-by-layer in-plane rumpling profile plotted in 
Fig. [3] reveals that the effect is highly localized at the p- 
type Ti02/K0 interface, quickly decaying upon moving 
into BaTiOs, and it is completely absent at the n-type 
BaO/Nb02 interfaces. 

The origin of the hidden interfacial in-plane polariza- 
tion can be easily understood with a simple electrostatic 
model based on formally charged ions. In Fig. [5l we 
compare the atomic structure of bulk BaTiOs [Fig.|5[a)] 
and KNbOs [Fig. [5|b)] with that present in the p-type 
Ti02/K0 [Fig. Ifc)] and n-type BaO/Nb02 [Fig. jsjd)] 
interfaces. The cleavage of bulk BaTiOs and KNbOs to 
form the p-type Ti02/K0 interface is accompanied by a 
change in the local electrostatic potential felt by the Ti 
cations. At the interface, some of the Ba cations (nom- 
inal charge +2) in the first neighbour atomic layer are 
replaced by K cations (nominal charge +1). This implies 
that at that interface the in-plane Ti cation movement is 
less constrained than in bulk, due to the reduced repul- 
sions with K ions with respect to Ba ones. Analogously, 
in the n-type interface the electrostatic potential felt by 
the Nb atoms is also altered. In this case the K+^ cations 
are replaced by Ba+^ ions, leading to an enhanced repul- 
sion that hinders the appearance of any in-plane polar- 
ization. As this effect is directly related to the nominal 
charges of the ions in both interfacial materials we would 
expect it to be present in other similar systems and could 
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be used to induce in-plane polarizations in other ferro- 
electric nanostructures. 



(a) 



+2 -2 +2 





+2 -2 +2 



(c) 



+2 -2 +2 





+1 -2 +1 +j -2 

Ba^^ ^ ^ O'- 

^ Nb^^ ^ 



FIG. 5. (Color online) Schematic view of bulk unit cells of 
(a) BaTiOs and (b) KNbOs, together with the ideal atomic 
structure at (c) p-type Ti02/K0, and (d) n-type BaO/Nb02 
interfaces. Atoms are represented by balls following the same 
conventions as in Fig. [3] Numbers indicate the nominal charge 
of the different ions. Solid blue (dashed red) arrows represent 
an enhancement (depletion) of the in-plane electrostatic re- 
pulsion between ions. The schema shows how the replacement 
of ions favors the movement of the transition metal B-cation 
in Ti02/K0 interfaces and hinders it in the BaO/Nb02 in- 
terfaces. 



To further characterize the development of the in-plane 
polarization, we take the relaxed coordinates of the asym- 
metric slab with m = 5 and / = 8 and remove by hand the 
in-plane displacements. We use this new structure as a 
reference configuration. Then, we compute the in-plane 
distortions required to go from the reference configura- 
tion to the relaxed structure and decompose into their 
X and y components. Finally, given fractions of these 
distortions are frozen in on top of the reference struc- 
ture. The energy surface obtained for different fractions 
of the in-plane rumplings is represented in Fig. [6j From 
this energy landscape we can extract (i) the stabiliza- 
tion energy (the energy difference between the polar state 



and the reference one), that amounts to 43.5 meV/slab, 
and (ii) the energy barriers that prevents the system to 
rotate the in-plane polarization from [110] direction to 
any of the other symmetry-equivalent positions, passing 
through the transition state at [100] positions. This is 
only 16 meV, slightly larger than those found in arti- 
ficial Ruddlesden-Popper-type superlattices.^^ To gauge 
the magnitude of these barriers, we compare them with 
the rotation zero-point-energy (ZPE) in this slab. This 
value is obtained through the nuclear motion Hamilto- 



H 



2M* 



dr]l 



dVy 



(12) 



where V (rjx^rjy) is the energy represented in Fig. [g] and 
M* can be shown to be 



1 



E 



(13) 



by writing the kinetic energy operator expressed in terms 
of atomic coordinates and masses (M^) using the linear 
trasformation relating the effective modes rjx and rjy with 
the atomic displacements through the coefficients q. To 
solve the Schrodinger equation associated to Eq. (IT2| and 
find the vibrational levels associated to Fig. [6] we follow 
the recipe given in Ref. 60, Then, we compute the ZPE 
by comparing the position of the first level to the mini- 
mum of the energy surface. The resulting value for the 
ZPE, 13.6 meV, is much smaller than the global stabi- 
lization energy but comparable to the height of the rota- 
tional barrier, indicating that quantum fluctuations will 
be important in these nanostructures. These fluctuations 
prevent the observation of any in-plane spontaneous po- 
larization, since the system will be delocalized over the 
four equivalent minima in [110] directions in a similar way 
to what happens in a dynamic Jahn- Teller problem.^^ 

This case can be compared with what happens in an 
incipient ferroelectric. In the latter case quantum fluctu- 
ations make the ZPE larger than the double well stabi- 
lization energy making the maximum of the probability 
density associated to the distortion (and the polarization) 
to be localized at the centrosymmetric state (origin in 
Fig.[6|. On the other hand, in the present case, the max- 
imum of the probability density corresponds with a non- 
zero polarization region around the origin (see dashed 
line in Fig. [6| but directional averaging results in a null 
net polarization. However, the coherent dynamics be- 
tween the wells could be disrupted by small electric fields 
along the plane, that would induce large changes in the 
directionality of the in-plane polarization. The signature 
of such potential energy surface for polarization rotation 
would be a high dielectric constant.^ 
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FIG. 6. (Color online) Two-dimensional energy surface as 
a function of the in-plane distortions of the BaTiOa atoms, 
as indicated in the main text. The dashed line follows the 
minimum of the valley. Units in meV. 



V. CONCLUSIONS 

In summary, using accurate first-principles simulations 
we have studied the influence of the ferroelectric polar- 
ization of the polar layer in the formation of 2DEG at 



BaTiOa/KNbOs interfaces. The most important conclu- 
sions that can be drawn are: (i) the spontaneous polariza- 
tion of the polar KNbOa layer cancels out almost exactly 
the "built-in" polarization discontinuity at the interface; 
(ii) as a consequence of this compensation, the critical 
thickness for the formation of the 2DEG is estimated to 
be between 42 and 44 unit cells of KNbOa, one order of 
magnitude larger than in SrTiOa/LaAlOa interfaces; (iii) 
this behavior can be easily explained in terms of a simple 
model based on the modern theory of polarization and 
basic electrostatics; and (iv) surprisingly, BaTiOa dis- 
plays an in-plane component of the polarization at the 
p-type Ti02/K0 interface, even when BaTiOs is under 
in-plane compressive strains. However we do not expect 
this in-plane polarization to be experimentally observable 
as the barriers for its rotation are very small and quan- 
tum fluctuations will prevent it from being localized in a 
particular direction. This situation could be easily mod- 
ified by the application of small electric fields that will 
break the symmetry and reveal the hidden polarization. 
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